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ABSTRACT 

The broad X-ray iron line, detected in many active galactic nuclei, is likely to be 
produced by fluorescence from the X-ray illuminated central parts of an accretion 
disc close to a supermassive black hole. The time-averaged shape of the line can be 
explained most naturally by a combination of special and general relativistic effects. 
Such line profiles contain information about the black hole spin and the accretion disc 
as well as the geometry of the emitting region and may help to test general relativity 
in the strong gravity regime. In this paper we embark on the computation of the 
temporal response of the line to the illuminating flux. Previous studies concentrated 
on the calculation of reverberation signatures from static sources illuminating the disc. 
In this paper we focus on the more physically justified case of fiares located above the 
accretion disc and corotating with it. We compute the time dependent iron line taking 
into account all general relativistic effects and show that its shape is of very complex 
nature, and also present light curves accompanying the iron line variability. We suggest 
that future X-ray satellites like XMM or Constellation-X may be capable of detecting 
features present in the computed reverberation maps. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - galaxies: 
Seyfert - X-rays: galaxies 
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The most extreme continuum variability in AGN is usually 
seen in the X-ray spectral band. The X-ray variability power 
spectrum reveals that low-amplitude fast flickering is an al- 
most permanent feature, but major changes, such as a flare 
with a large change in luminosity (like doubling) occur much 
less frequently. The duration of such flares may by of the or- 
der of a few thousand seconds (e.g. see figures in McHardy 
(1989)). Recent spectroscopic observations suggest that such 
X-ray flux variations can be accompanied by Fe Ka line vari- 
ability which occurs on a similar time scale (Iwasawa et al. 
(1996)). This might indicate that the size of the X-ray con- 
tinuum emitting region is small in comparison with other 
sources contributing to the overall continuum and that the 
iron line is produced in a compact region close to the con- 
tinuum source and is correlated with the X-ray outburst. 
The X-ray emission is thought to originate from the inner- 
most part of the accretion disc. This region is subject to a 
number of instabilities which may produce flare- like activity. 
There is a large body of observational and theoretical work 
which supports this view. Among the proposed scenarios 
invoked to explain such a phenomenon are: the ampliflca- 
tion of magnetic fields in the disc by convection and differ- 



ential rotation leading to the emergence of magnetic flelds 
carrying hot plasma (Galeev, Rosner & Vaiana, 1979), hot 
magnetic arcs due to the Parker instability (Chagelishvili, 
Leminadze & Rogava (1989)), magnetic flares due to explo- 
sive release of stored magnetic energy in the accretion disc 
(Vries & Kuijpers (1992)) and vortices (Abramowicz et al. 
(1992)). The fluorescent iron line emission may be produced 
by irradiation of the disc material by such flares located 
above the accretion disc (Lightman & White (1988), George 
& Fabian (1991), Matt, Perola & Piro, (1991)). As a resufi, 
the observed time-averaged spectra have distinctive, skewed, 
double-peaked profiles which reflect the Doppler and grav- 
itational shifts in a strongly curved spacetime (Fabian et 
al. 1989, Laor 1991, Mushotzky et al. (1995), Tanaka et al. 
(1995), Nandra et al. (1997)). Apart from the line emission, 
the disc material should also reflect a part of the X-ray con- 
tinuum. Evidence for this effect was found by Pounds et al. 
(1990) for a number of Seyfert galaxies. Since the spectral 
properties of the Fe Ka line are well known, the iron line 
may be used as a probe of the very near environment of a 
black hole and the black hole itself. This may enable one to 
study the geometry of the emission region, estimate black 
hole masses or search for the dragging of inertial frames. 
The frame dragging effect is especially interesting because 
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its presence, in cooperation with the shear viscosity of an 
accreting material, may be responsible for the alignment of 
an accretion disc symmetry axis with the spin axis of a black 
hole or help to explain the stability of jets (Bardeen & Pe- 
terson (1975)). 

Much theoretical work has been done to search for the 
signatures of the mass and spin parameter of black holes 
as well as the geometry of the line and continuum emitting 
regions. Laor (1991) realized the importance of the impact 
of general relativistic effects in the vicinity of the rotating 
black hole on the line shape. This influence of the black hole 
spin on the iron line was later quantified by Dabrowski et 
al. (1997). They calculated the iron line profile assuming 
that the line is produced in the thin corona above the accre- 
tion disc. They considered the case of the Seyfert 1 galaxy 
MCG-6-30-15 and concluded that the very broad and skewed 
iron line can be explained if the emission originates from 
the disc extending down to the radius of marginal stabil- 
ity Tms < 1.23m which corresponds to the spin parameter 
a > 0.94. Thus they claimed that the black hole is rapidly 
rotating and should be described by Kerr geometry. A simi- 
lar conclusion was drawn by Bromley, Chen, & Miller (1997) 
who suggested that the iron line profile was more likely to 
originate from a disc about a rotating black hole. Another 
approach was presented by Reynolds & Begelman (1997) 
(hereafter RB97) who were able to explain the same data 
by assuming Schwarzschild geometry (i.e non-rotating black 
hole) and including the matter below the innermost stable 
orbit. However in their model the iron line was excited by an 
external flare located on the symmetry axis above the accre- 
tion disc and illuminating the regions below the innermost 
stable orbit. They also claimed that the rapidly rotating 
black hole in the radio-quiet active galaxy MCG-6-30-15 is 
contrary to some theoretical models invoked to explain radio 
loudness in AGN (e.g. Rees et al. (1982), Wilson & Colbert 
(1995)). However more recent studies show that this is not 
necessarily the case (Ghosh & Abramowicz (1997), Livio, 
Ogilvie & Pringle (1999)). Recently Young, Ross & Fabian 
(1998) calculated the spectrum corresponding to the RB97 
model and found that it should possess a photoelectric ab- 
sorption edge of iron, which is not seen in the data. However 
with the quality of the current data the issue of the black 
hole spin remains ambiguous. An entirely different approach 
to searching for the observational signatures of the black hole 
spin, mass and the geometry of emission region was first sug- 
gested by Fabian et al. (1989) and later considered in any 
detail by Stella (1990). He computed the iron line variabil- 
ity from a disc around a Schwarzschild black hole assuming 
a pointlike variable X-ray source located in the geometric 
centre of the disc. For an observer at infinity any such vari- 
ations of the primary X-ray source would be 'echoed' by 
different locations on the disc leading to time evolution of 
the line profile. The primary goal of this research was to 
suggest a method to estimate black hole masses. His work 
was extended by Matt & Perola (1992) and Campana & 
Stella (1993, 1995) who considered modified source geome- 
tries. Recently Reynolds et al. (1999) (hereafter RYBF99) 
generalized these results by assuming arbitrary black hole 
spin and searched for observational signatures of this pa- 
rameter. They mainly focused on the case of on-axis source 
but briefly discussed preliminary results for a static flash- 
like off-axis flare in the vicinity of a Kerr black hole. 



As noted above, the duration of bright flares can 
be of the order of lOOmr in geometric time units 
(IGmrlO^Mo/c^ « mrAQs; for the black hole mass M = 
lO^m-rMo). It is very likely that any variability associated 
with blobs above accretion discs most likely comes from 
moving sources. Thus in general, the assumptions of a static 
primary X-ray source made in previous studies, and in par- 
ticular the assumption of an (5-like flash from a static flare 
made by RYBF99, limit the applicability of their results. In 
the present work we relax these assumptions and consider 
reverberation effects from corotating flares above an accre- 
tion disc. We take into account the effect of Doppler boost- 
ing of radiation from a moving flare which was neglected 
by RYBF99 and allow the flares to revolve around a black 
hole. We calculate the time dependent iron line for different 
values of inclination of the accretion disc and spin parame- 
ter (we consider the Schwarzschild and the almost maximal 
Kerr cases) and different positions of the flare relative to the 
accretion disc. We also compute the light curves accompa- 
nying the iron line variations including the direct flux from 
the flare and the reflected flux from the disc. Similar calcu- 
lations were performed for example by Bao (1992), Karas et 
al. (1992), Zakharov (1994) and Bromley et al. (1997) in the 
context of corotating spots located on the accretion disc. 

The following Section and its subsections describe the 
assumptions made in our model and the algorithm used to 
simulate the light curves and the iron line flux variations. 
The last Section is devoted to the presentation and discus- 
sion of our results. 



2 THE MODEL ASSUMPTIONS AND THE 
ALGORITHM 

2.1 The X-ray flare and the disc illumination 

We assume that the primary source of radiation is an X- 
ray flare situated above the accretion disc. As noted in the 
introduction, magnetic instabilities may cause such a flare 
to emerge from the disc and its position relative to the disc 
surface may remain approximately unchanged for some time. 
Thus we make the assumption that the flare is confined by 
the magnetic field and corotates with the disc and possesses 
fixed Boyer-Lindquist coordinates r and 0. We also assume 
that: the flaring region is point-like, has a finite life-time and 
that photons propagate freely to either the observer, to the 
black hole, to the disc or escape to infinity once emitted by 
the blob. This means that the corona above the accretion 
disc is optically thin. This is a generalization of the con- 
siderations presented by RYBF99 who model X-ray flare as 
a (5— like impulse and additionally assume that radiation is 
produced in the locally non-rotating frame of reference. In 
order to numerically integrate the photon trajectory from 
the corotating flare we first analytically derive appropriate 
constants of motion to propagate photons through the Kerr 
metric. All the necessary formulae are given in Appendix 
B, where we express the constants of motion in terms of 
polar and azimuthal angles in the local rest frame of the 
point-like source. Thus we can easily model the isotropic 
distribution of radiation in the emitter's frame by a Monte 
Carlo method. This enables us to calculate the illuminating 
flux as a function of time and energy in the rest frame of the 
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corotating disc material taking into account all the general 
relativistic effects, including previously neglected Doppler 
boosting from the moving flare. In the implementation of 
the algorithm wo used the formula (see Appendix C): 



F{Ei,isd)<^SE^=gl 






fsd 

'yd<j}dr 



1/2 



where we additionally assume that the flare emits a power 
law spectrum with energy index a. The factor fsd, which 
is the ratio of the number of photons intersecting a small 
patch on the disc (defined by d(jf> and dr) to the total number 
of emitted photons, was calculated by means of the Monte 
Carlo method. In the above formula tsd is the time it takes 
for photon emitted from the source to reach the disc element 
and 7 is the Lorentz factor for the relative motion of the 
disc element and the locally non-rotating observer. Qsd is 
the redshift factor given by 



9sd = 






where pdp and ps^ are the photon four-momenta at the disc 
and source respectively. We follow the lines of reasoning of 
Cunningham (1975) and split the velocity field of the accre- 
tion disc Ud into the region outside the radius of marginal 
stability, where matter follows circular orbits, and the re- 
gion within the innermost stable orbit where the matter has 
a negative radial component and spirals towards the black 
hole. The formulae for the flare velocity field u^ are collected 
in Appendix B. 



2.2 Reprocessing of radiation by the accretion 
disc 

The radiation which impinges upon the disc surface is re- 
processed and reflected. The shape of the reflected spectrum 
depends upon many factors, the most important being the 
disc structure and ionization state of the accretion disc. In 
the subsequent discussion we assume the disc to be opti- 
cally thick and geometrically thin. The optical depth of the 
accretion disc was calculated by many authors (e.g. Frank, 
King, & Raine. (1995), RB97) and is shown to be large. This 
guarantees that the illuminating radiation can be efficiently 
reflected. However the albedo of the reflecting material de- 
pends on the ionization parameter which varies across the 
disc surface: (,{r,(j)) = 4n Fx{r,(f))/[n{r) cos n], where n{r) 
is the comoving electron number density given by Eq. 4 of 
RYBF99, cosn is the angle between the incident rays and 
the normal to the disc surface in the comoving frame of the 
disc and Fx is the X-ray flux deflned over some flxed energy 
band at a particular position on the disc. We calculate the 
parameter ^ and use the results of Ross, Fabian & Young 
(1998) to calculate the albedo aaiE) as a function of energy 
for all grid elements on the disc. We then adopt the simple 
prescription of RYBF99 to separately model the iron line 
fluorescence in the rest frame of the moving disc element, 
i.e. we assume that for ^ < 100 erg cm s^^ the matter emits 
cold iron line at 6.4 keV with some yield Y , for ^ in the 
range from 100 erg cm s^^ to 500 erg cm s~^ there is no 
emission because the radiation is destroyed by Auger mech- 
anism, for ^ between 500 erg cm s~^ and 5000 erg cm s~^ 
we use a composition of two lines at 6.67 keV and 6.97 keV 
each with the fluorescent yield Y and assume that no iron 



line is produced for ^ > 5000 erg cm s^^ because of total 
ionization of the matter (here we approximate the iron lines 
as delta functions) . Note that these additional complications 
are essentially relevant only for the Schwardschild black hole 
because only in this case can a significant fraction of the disc 
material within the innermost stable orbit (vms = Gm) be in 
the higher ionization state (large 5). This arises from the fact 
that the radial density gradient of the matter spiralling to- 
wards the black hole is strongly negative in this region so the 
disc significantly decreases its density whereas the density 
changes for radii greater than rms are comparatively negli- 
gible and the disc remains relatively dense and thin, at least 
for small accretion rates. This means that the disc may be 
treated as approximately cold for r > Vms and the variations 
of the ionization parameter induced by the density changes 
should not dramatically alter the reflection from the disc in 
this region. Thus the assumption of RYBF99 about the high 
density and negligible ionization of the region outside Vms in 
the Schwarzschild case may be justified. In the Kerr case the 
radius of marginal stability asymptotically approaches the 
event horizon and we assume that the disc may be approx- 
imated as being cold on average. The results of Rozaiiska 
et al. (1998), who show that the presence of a hot corona 
makes the disc denser and less ionized, further strengthen 
the above argument. 

In the subsequent considerations we use the additional 
simplifying assumption of isotropic refiection and write the 
reflected intensity in the form: 

lR{Ed,isd) = -aR{^,Ea)F{Ea) 



2.3 Calculation of light curves and iron line 
variability 

In order to calculate the time-varying iron line profiles and 
the corresponding continuum light curves we used the fol- 
lowing prescription: 

1) We use a ray-tracing technique to follow the trajecto- 
ries of photons from positions (x, y) on the observer's image 
plane to the disc, keeping track of all Boyer-Lindquist coor- 
dinates idoix,y), r{x,y), 6{x,y), (l){x,y), until the photons 
either intersect the accretion disc plane or disappear below 
the event horizon (The observer is located at r = 1000m, 
6 — i and = 0°, where i is the inclination) . We then calcu- 
late the redshift factor corresponding to the particular final 
position of a photon and its arrival time from the disc. We 
generate an image of approximately 1600 x 1600 pixels cov- 
ering roughly an area of 1607n x 160m. 

2) Using the time dependent ionization and illumination 
pattern described above and albedo aii{E, ^), we specify the 
emergent fiux from the disc as a function of time and en- 
ergy in the fixed spectral band, i.e. we integrate all photons 
coming from regions of equal arrival time regions on the ac- 
cretion disc at a particular energy. In doing so we also take 
into account the fact that the source changes its position 
and rotate the illumination and ionization patterns accord- 
ingly. This allowed us to compute the reflected spectrum as 
a function of time and thus the reflected component of the 
light curve. We use the following formulae (see appendix C 
for additional explanations): 
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Fn{t)^J2^n\i), 



where Fj^ {t) is the contribution to the flux from the flare 
at the n-th position given by; 



F'-\t) 
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where Qdo is the redshift factor for the disc-to-observer 
case, dxdy is the pixel surface area on the observer's im- 
age plane, ro is the distance from the black hole system and 

3) We compute the direct flux component by propa- 
gating photons from the source to the observer. We use a 
formula analogous to that above to calculate this flux. 

4) Having specified the illumination pattern and ioniza- 
tion state for a given position of the flare above the accretion 
disc, we determine the contribution of each grid element on 
the disc surface to the iron line profile from the formula: 

1 



*("'(£,i)(X 



AtAE 



E 



Y{^,E/gdo)gLFx{r, 



t^t+At,E^E+AE 



where Fx (r, (f)) is the integrated over a fixed energy band X- 
ray fiux at the given position on the disc. The above function 
is the Green's function (transfer function) corresponding to 
the n-th position of the flare. 

5) In the last step we use the transfer functions for 
single isolated flares in order to calculate the observed iron 
line flux variations. We use the following expression: 



Ffe(£,t)0C^*W(i5,t), 



where we also take into account the delays due to the rev- 
olution of the flare around the black hole. As in the case of 
the reflected flux, we interpolate the contributions to the flux 
from the isolated flares and then add the interpolated fluxes. 
Note that the binning of the transfer functions in time would 
introduce a small additional (artiflcial) flux variability, be- 
cause we physically rotate the illumination pattern. This is 
because the differences between the arrival times from the 
neighbouring disc elements on the approaching side of the 
disc would be shorter compared to the corresponding differ- 
ences on the opposite side. 



3 RESULTS 

We present iron line flux variations and the corresponding 
light curves for one full revolution of the flare around the 
black hole with a step-like luminosity variation in the local 
frame of reference of the flare. This is for illustrative pur- 
poses and we stress that while it is trivial, for example, to 
generalize the code to include arbitrary intrinsic luminosity 
variations in the rest frame of the flare, we do not wish to 
complicate the interpretation of the results. Future observa- 
tions may guide us as to how to introduce additional modifl- 
cations in the computations. The results in the present form 
can in principle be applied to the case of isolated flares or 
at least strong outbursts accompanied by minor ones. In the 
computations we assume that the source emits a power law 
spectrum with energy index a = 1. 
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Figure 1. Total light curve (solid line), the direct flux from the 
flare (dot-dashed line) and the small contribution from the re- 
flected component (dashed line) for the same set of parameters 
as in Fig. 3. The flare is switched on from = 0° to (^ = 360°. 
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Figure 2. Total light curve (solid line), the direct flux from the 
flare (dot-dashed line) and the small contribution from the re- 
flected component (dashed line) for the same set of parameters 
as in Fig. 4. The flare is switched on from = 0° to (^ = 360°. 



Fig.l and Fig. 2 present the light curves. The light curves 
are a superposition of the direct flux from the flare and the 
reflected flux from the disc. The small offset between these 
two components (note the shift between initial and final 
times of the reflected and direct flux) is due to the addi- 
tional time photons need to propagate towards the disc and 
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Figure 3. The colour map represents time-sequence of tlie iron 
line spectrum, i.e. log(F(time, energy)), where F is the flux di- 
vided by the maximum flux during the whole time interval. 
The Boyer-Lindquist coordinates of the flare are: Ofiare = 70°, 
^ flare = 10m with 4' flare, changing from 0° to 360°. The spin 
parameter and the inclination are a = 0.998 and i = 30° respec- 
tively. 



Figure 4. As Fig. 3, but for i = 50°. 
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back. The small albedo of the disc material in the chosen en- 
ergy band (2-9 keV) results in a relatively small contribution 
from reflection continuum to the total light curve. However 
the albedo strongly increases with energy and thus the reflec- 
tion would be much more efficient in the harder part of the 
observed spectrum. Another effect which would enhance the 
reflected component (and suppress the direct one) is the fo- 
cusing of light from the flare. More radiation from the flares 
closer to the central black hole would be focused towards the 
disc and reflected, and simultaneously less direct flux would 
reach the observer (Martocchia & Matt, (1996)). The peaks 
in the direct light curve are mainly attributed to the Doppler 
effect. Note how the flux variability is enhanced by the in- 
clination of the system. For higher inclinations the Doppler 
effects are more pronounced and the peaks are narrower and 
have higher maxima. It is important to note that the flux 
variability observed in the chosen spectral band (2-9 keV 
in our case) is also a function of the energy index and will 
be increased in the case of the steeper spectra. The peaks 
due to gravitational focusing (e.g. Cunningham & Bardeen, 
(1973)) are not seen here because we limit our discussion to 
low inclinations in order to avoid additional complications 
related to the limb darkening. This is not a serious limitation 
to our model because the majority of Seyfert 1 galaxies are 
likely to have small inclination angles. However we compare 
our results for the direct flux in the test case of a source 
orbiting the black hole in the equatorial plane for very high 
inclinations of the orbit and find agreement with the results 
of Cunningham & Bardeen (1973). 

Figs from 3 to 8 show the evolution of the iron 
line spectra after the subtraction of the varying power law 




Time (Gm/c') 



Figure 5. As Fig. 3, but for a = 0.0. 

continuum for different sets of parameters. As for the case 
of the light curves, the time offset of the initial spectrum is 
the difference between the arrival of the direct contirmum 
flux and the first appearance of the iron line. There is a 
clear difference between the results for the two considered 
extremal values of the spin parameter (e.g. compare Fig. 4 
and Fig. 6). The iron line extends to lower frequencies in the 
Kerr case which is due to the cold iron line being produced 
very close to the black hole. The influence of the ionization 
of the accretion disc below the radius of marginal stability 
rms = 1.23m. is relatively unimportant. On the contrary, in 
the case of the Schwarzschild black hole the bulk of the cold 
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Figure 6. As Fig. 3, but for a = 0.0 and i = 50°. 
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Figure 8. As Fig. 3, but for rfiare = 3m. 
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Figure 7. As Fig. 3, but for rfiam = 6.5m. 

iron line is produced only above Vms ~ 6m and accretion disc 
may be in highly ionized within this radius. Thus for r < 6m 
the line may be either destroyed by the Auger process or not 
produced because of the total ionization of the disc material. 
The radiation from this region would be subject to high en- 
ergy shifts which means that the lack of the iron line flux at 
lower energies (see Fig. 5 and Fig. 6) is the result of the fact 
that some parts of the disc below rms = &rn cannot produce 
the iron line. For the flare at the position (fifiare the regions 
on the disc responsible for the production of the iron line in 
the higher ionization states (500 erg cm s~^ < C < 5000 erg 
cm s~^) are located around (j> — 0/iare±9O°. The regions on 
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Energy (keV) 

Figure 9. Integrated iron line profile for t = (70 H- 80)Gm/c'^ 
after the detection of the first signal from the flare for: a = 0.998, 
9 flare = 70°, r flare = 10m and i = 30° 



the opposite side of the black hole (0 — (j) flare + 180°) and 
nearest to the flare {(j> = 4' flare) have even higher ionization 
parameter which is the result of the strong light bending 
and the direct illumination of the disc surface respectively. 
This effect gives rise to the varying width of the line for 
a = with a minimum around (100 -^ 150)Gm/c^ (see Fig. 5 
and Fig. 6). Future high-throughput instruments like XMM 
and Constellation-X may be capable of detecting such fea- 
ture in the temporal iron line profiles. However we stress 
that a detection of this weak feature would be difflcult be- 
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Figure 10. Integrated iron line profile for t = (70 -H 80)Gm/c'* 
after the detection of the first signal from the fiare for: a = 0.998, 
d flare = 70°, rjiarf, = 10m and i = 50° 
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Figure 12. Integrated iron line profile at i = (50 -H 55)Gm/c'^ 
after the detection of the first signal from the fiare for: a = 0.998, 
& flare = 70°, r flare = 3m and i = 30° 
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Figure 11. Integrated iron line profile at i = (40 -^ 45)Gm/c^ 
after the detection of the first signal from the flare for: a = 0.998, 
& flare = 70°, r flare = 6.5m and i = 30° 



cause it could be overwhelmed by the continuum variations. 
Other interesting features are the two 'bumps' seen for the 
inclination i = 50° and the Schwarzschild black hole. Here 
the iron line extends to higher frequencies than the max- 
imum frequency in the Kerr case for the same inclination 
of the disc (compare Fig. 4 and Fig. 6). These features are 
due to the highly blueshifted 'hot' iron lines {E\.J^^ = 6.67 



keV and E. 



(2) 



6.97 keV) from the ionized regions. This 



effect should also be within the reach of capabilities of XMM 
and particularly Constellation-X. It has to be stressed how- 
ever that the above results are sensitive to the parameters 
and depend on the position of the flare and its X-ray ef- 
ficiency (here assumed to be r]x = 0.06) and the density 
of the disc below the innermost stable orbit (which is regu- 
lated by the half thickness of the disc hdisc', here we assumed 
hdisc = 0.05r). In principle, a source of very low efficiency 
located at a very large distance and large 6 fiare and illumi- 
nating the high density region below the innermost stable 
orbit could not ionize it, and a faint cold iron line could also 
be produced all the way down to the event horizon in the 
Schwarzschild case. In such situation it would be difficult to 
distinguish between the Kerr and Schwarzschild black holes. 
On the other hand, the flares located closer to the black 
hole and closer and less inclined with respect to the rota- 
tion axis would ionize the regions below the innermost orbit 
much more efficiently leading to even more pronounced dif- 
ferences (e.g. then those shown in Fig. 4 and Fig. 6) between 
the temporal iron line profiles for the two extreme values 
of the spin parameter. The work of RB97 and Young et al. 
(1998) has shown that in the special case of a nonrotating 
black hole and the flare located on the rotation axis very 
close to the black hole (e.g. h = 3.5m) essentially all radia- 
tion is focused towards the disc (and the black hole) leading 
to strong ionization. Young et al. (1998) also demonstrated 
that in this case the iron absorption edge becomes impor- 
tant. This effect leads to the trough in the iron line spectrum 
above the considered iron line rest energies. Thus more ac- 
curate modelling should include the absorption edge which 
would decrease the iron line flux from the ionized regions be- 
low Tms and create an oscillating absorption feature in the 
reverberation maps. It is also possible that other processes 
could influence the iron line production below rms- For ex- 
ample, the radiation drag could decrease the momentum of 
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the infalling matter and thus accelerate the accretion be- 
low Vms leading to a decrease in density and an increase in 
the ionization parameter in this region. Thus, as a result of 
the above processes, the iron line in the Schwarzschild case 
may not extend to frequencies as low as those in the Kerr 
case and may possess the absorption features, so in general 
the appearance of the reverberation maps for the two cases 
should indeed be different, making it possible to distinguish 
between a rotating and nonrotating black hole. 

The most prominent feature seen in all diagrams show- 
ing the evolution of the iron line spectrum is the drifting 
maximum of the flux. This sinusoid-shaped feature is due to 
the fluorescence from the most strongly illuminated part of 
the disc located just below the flare. Its shape and bright- 
ness is the result of special and general relativistic effects. 
Initially the flare moves in a direction perpendicular to the 
line of sight and the maximum is redshifted due to the trans- 
verse Doppler effect and gravitational frequency shift. Later 
it gradually begins to recede from the observer. The bright- 
ness of this maximum is suppressed at this stage because the 
flux from the highly illuminated part of the disc below the 
flare decreases as a result of light beaming in the opposite 
direction relative to the observer. The other interesting fea- 
ture seen in this phase is the double maximum clearly seen 
around (50 ~ 100)Gm/c'^ in Fig. 4 and Fig. 6. As mentioned 
above, the lower maximum is the result of the redshifting of 
the bulk of the radiation from below the receding flare. The 
emerging upper maximum on the other hand, is due to flu- 
orescence from the illuminated regions on the approaching 
side of the disc. As the flare rotates and enters the approach- 
ing side of the disc the bulk of the iron line flux shifts towards 
high energies and the Doppler boosting of the radiation re- 
flected by the strongly illuminated regions of the disc leads 
to the rise of the iron line maximum and its brightening. Af- 
ter the flare passes through cpfiare ~ 270" the energy shift 
decreases and gravitational and transverse Doppler shifts 
gradually take over. Note that during the whole revolution 
of the flare the redshift factor g < 1 for a longer time, be- 
cause the gravitational and transverse Doppler shifts may 
for some time overwhelm the blueshifts when the flare is on 
the approaching side of the disc, i.e. when the flare has just 
entered or is just about to leave the approaching side of the 
disc. This effect gives rise to a flattening of the minimum (in 
energy) of the main flux maximum, i.e. the bright A-shaped 
feature appears to be sharper (e.g. see Fig. 4 and Fig. 6). 

A very interesting feature can be seen when the flare 
is receding {(f> « 90°). The main redshifted flux maximum 
broadens (see Fig. 4 and Fig. 6) or even splits into two max- 
ima (see Fig. 3, Fig. 5 or Fig. 7). In order to illustrate this 
effect better, appropriate cross-sections through the rever- 
beration maps shown in Fig. 3 and Fig. 4 are given (see 
Fig. 9 and Fig. 10 respectively). Note that in the case of 
Fig. 9 three maxima are seen. This thickening or splitting 
is a purely general relativistic effect. The middle maximum 
is a consequence of strong bending of light which is focused 
on the opposite side of the black hole relative to the actual 
position of the flare and leads to the enhancement of the 
illumination of the disc in this region. The influence of this 
effect is even more important for the flare located closer to 
the black hole {vfiare = 6.5m) as shown in Fig. 7 and Fig. 
11. Fig. 11 shows three distinct maxima with the middle 
maximum due to light focusing being dominant. As in the 



previous cases, such an effect could be observable by future 
high throughput spectrometers. The magnitude of the en- 
ergy shift of the main maximum is related to the inclination 
of the accretion disc (e.g. compare Fig. 3 and Fig. 4) and also 
to the distance of the flare from the centre. The latter effect 
is demonstrated for example in Fig. 3, Fig. 7 and Fig. 8 where 
the bright A-shaped feature extends gradually to lower en- 
ergies as the distance rfiare of the flare from the centre de- 
creases. The same figures also illustrate how the timescale of 
the sinusoidal variation changes as a function of rfiare with 
shorter time scales corresponding to faster orbital motion of 
the closer flares. Note how the structures on the reverber- 
ation map change for a very close flare {rfiare ~ 3m; see 
Fig. 8). There is no triple maximum on the left hand side 
of the A-shaped feature, but the maximum due to focusing 
emerges on the opposite side. This is also shown on Fig. 12 
where three maxima are visible but, unlike in the case of 
more distant flares, in this case the lowermost and middle 
maxima are due to gravitational focusing. The lack of the 
uppermost and middle maxima on the left hand side of this 
feature is related to the very large overall gravitational and 
transverse Doppler redshift. It is also caused by the faster 
motion of the flare (and thus the main flux maximum) com- 
bined with the relatively slow propagation of signals slowed 
by the Shapiro delay in the vicinity of the black hole. In 
other words, the ratio between the time it takes for the flare 
to complete a half of one full revolution around the black 
hole to the signal crossing time decreases with decreasing 
distance of the flare from the centre. The above effects con- 
tribute to the swamping of these two maxima by the main 
bright maximum. 

We have demonstrated that the very complex behaviour 
of the lines is predicted by a relatively simple model of the 
orbiting flare. Our results can be applied to future obser- 
vations made with XMM and particularly Constellation-X 
and may help to understand various processes operating in 
the vicinity of black holes. These results can be used to put 
constraints on their masses and spin parameters, the geom- 
etry of the emitting region and test general relativity in the 
strong gravity regime. Particularly interesting features seen 
on our diagrams are the triple maxima in the temporary 
iron line profiles. The middle maximum is a consequence 
of strong light bending by the extreme gravitational fields 
very close to the central black hole and is thus a prediction of 
general relativity. This may indicate that a complicated tem- 
poral behaviour of the iron line profiles seen in future data 
may not necessary imply a complicated model but could be 
explained for example in the framework of the orbiting flare 
model. 
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APPENDIX A: THE KERR METRIC 

Space-time in the vicinity of a rotating black hole is de- 
scribed by the Kerr metric. In Boyer-Lindquist coordinates 
this metric reads: 



as = — I 1 — — I at ^ sm 

-f— sm 6a(p , 
where 



A 



A = r^ + a^ - 2r 



A = {r + a ) — a A sin 



APPENDIX B: THE CONSTANTS OF MOTION 
FOR THE OFF-AXIS SOURCE ON A 
SPATIALLY CIRCULAR ORBIT 

The photon trajectory may be specified by two constants 
of motion (the component of angular momentum parallel 
to the symmetry axis I and the Carter constant Q) which 
can be expressed in terms of the direction cosines ej of the 
photon momentum k and the comoving tetrad Aj. 

k ■ A; . „ 



(k-h-k)i/2 



where g is the Kerr metric tensor and h is the transverse 
projecting operator defined by h = g + u ■ u and u denotes 
the four velocity of the source given by: 



u = c{dt + na^, 

where 

2r 



C 



1 — ^(1 — afisin 6) — (r + a ) Q sin ( 



-1/2 



and Q is the angular velocity of the fiare deGned by: 



a +(r sin 61)3/2 



Note that such a definition does not violate (for the param- 
eters considered in Section 3) the obvious condition that the 
source must follow a timelike worldline: 



(u ■ u 



[gtt + "2^94,1 + ^ g,p4,] < 



The necessary components of the comoving tetrad and pho- 
ton momentum are: 

A,= (|) dr 

and 



fc"- = ±i [(2r - alf + A{r{r + 2) - L)] ^'^ 

1/2 



"E 

4u 



I'' 



2 2 
+ a cos ( 



A 

AE 



sin 
/E-2r 



1/2 



where u> = 2ar/A is the angular velocity of the frame drag- 
ging and L = Q+l^. The directional cosines can be expressed 
as functions of polar '^ and azimuthal $ angles in the rest 
frame of the source in the usual way: 

ef = cos ^ 



T, — r + a cos ( 



eg — sin \l/ cos $ 
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Here we chose the local z axis to point in the dr direction 
and the local x axis in de direction. We can now combine the 
above equations to obtain the desired set of two equations 
for the two constants of motion I and Q: 



S(e? + e^)k-h-k+ ^ 



sin^ e 



{2r -alf 
A 



+ r{r + 2) (Bl) 



2 
Sf , 2 2 , 

-rria cos ( 
ei 



-^^L)+L^ ^ ^ + r(:r + 2) (B2) 

in d A 



All components containing L in Eq. Bl cancel out mirac- 
ulously and we obtain relatively simple quadratic equation 
for I which has the solution: 



^ 1 / .^N 



X) \X 



1/2 



with 



<iv{^) ( ^ sin^ < 



+ C^[(fi-a;)g + p(f7)a;]- 



where dile = 4:n{dNEc/Ntot) is the solid angle in the rest 
frame of the emitter which corresponds to the surface ele- 
ment dS on the disc, dNse is the number of photons emit- 
ted into this solid angle, gsd is the redshift factor and Le^ 
is the luminosity of the source per unit energy. Note that 
in special relativity the transformation law for the solid 
angle is dOe = g^^dftobs- This means that in this regime 
we can recover the usual result Ftot oc gj^ by integrating 
the formula for the flux received by the disc element. Since 
the surface density of photons is locally constant we have 
dS = {dNEjdNEjdS, where dS = ■y{A/Ay'''^d(t>dr is the 
surface area of the fixed grid element on the disc and 7 is 
the Lorentz factor of the relative motion of the disc element 
and the locally non-rotating observer. For r < rms, where 
Vms is the radius of marginal stability, this factor is given 
by: 

^^(l_yW)2-,-l/2^ 

where 



ne = l/{r^^^ + a), 



1 fa 



A 



/(*,<&) 



Y = 2(-^] [lop{lo) + C'^ { ^Lo{Q - Lofq sin^ e - p\Q,)lo 



vae; 



4ar 



^ = "^ (^ + '■^^ + ^^ + "' °°'' ^) ^^*' "^^ 

+ (^)%H(cVc.)-i), 

where 

2r A . 2„ ./T ^N sin^vl/sin^^ 

p{x) = l- — + -u;xsm^e, /'!',$=- . 3, . 2^ 

1j 2-, 1 — sm W sm <P 

and g = 1 — (2r/E). The second constant of motion Q — 
L — I can be easily calculated from Eq. B2. 



APPENDIX C: ILLUMINATION OF THE DISC 
AND LIGHT CURVE FORMULAE 

The observed flux consists of two elements: the direct flux 
from the flare and the reflected flux from the accretion disc. 
Let us first consider the observed flux F{Ed) as seen by the 
disc element moving with the accretion flow. 



F{Ei) = 



EddNE^ 
dS '' 



(CI) 



where Ed is the photon energy in the matter rest frame and 
dNE^ is the number of photons per erg per second which 
illuminate an element of surface dS on the accretion disc. 
Because the number of photons has to be conserved and 
dtddEd ~ dtedEe we have dNE^ = dNE,, where the sub- 
script e and d denote the rest frame of the flare and the disc 
frame respectively. Thus Eq. CI takes the form: 



F{Ed) = 



gsdLE^ rffie 

dS 1^ 



whereas for rn < r < rms, where th is the radius of the 
event horizon: 

^^(l_yM2_^(«2-,-l/2^ 

where 



V 



(r) 



A u* 



y'-^^l'--. 



The variables u*, u^ and u'^ are the components of the 
4-velocity below the innermost stable orbit (see e.g. Cun- 
ningham 1975) and V^ , V^ and V'^ are the components of 
3-velocity relative to the locally non-rotating frame which 
can be calculated from V ^' = {u^e]i )/{u'^e]^ ), where 
j — r,4> and ej^ are the components of the basis of one- 
forms corresponding to the set of basis vectors of the lo- 
cally non-rotating observer (see Bardeen, Press, Teukolsky, 
1972). Thus for the luminosity in the form of the power law 
Lec oc E^" we get: 



F{Ed) ex Se, = g'J-E^' 



fsd /A 



■ydcpdr 



(^) 



1/2 



where fsd = dNE^/Ntot is the fraction of Ntot which cor- 
responds to the grid element defined by d<j} and dr. Assum- 
ing isotropic reflection, the final expressions for the reflected 
component of the flux read: 

n 

where Fj^{t) is the contribution to the flux from the flare 
at n-th position given by: 



^^"'® = if E C!{EouE^2,a) 



gdolR{Eo/gdo)dEo 



dxdy 



where t is the sum of the time tsd it takes for the photon 
emitted at the source to reach a given element on the disc 
and the photon arrival time tdo from that particular point on 
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the disc to the observer. Qdo is the redshift factor for the disc- 
to-observer case, C is the proportionahty coefficient which 
depends on the energy index and the observed spectral band, 
At is the width of the time bin and 

lR{Eo/gdo) = -aR{^,Ei)F{Ed), 

TV 

where a_R(^,i5) is the albedo of the reffecting material. We 
used the transformation law for the intensity Ieo = gtd^Ec 

(n) 

in the expression for F^ above. Note that in the formula 
for -FJj" we include contributions from all photon isodelay 
regions on the disc. We then interpolate Fj^ at a given 
time t for all flare positions and sum these contributions to 
get the reflected flux FR{t). We use a very similar method 
to calculate the iron line spectra (appropriate formulae are 
given in the main text). 

Analogously, we obtain the formula for the direct com- 
ponent of the flux in the following form: 

where fso is the fraction of the total number of emitted 
photons intersecting the observer's image plane which has 
the surface area Ax/S.y. Note that the formulae for Fr, and 
Fd have the same proportionality coefBcient which enables 
us to compare the relative flux contributions. 
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